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Abstract 

In this paper we propose a new numerical method for solving stochastic differ- 
ential equations (SDEs). As an application of this method we propose an explicit 
numerical scheme for a super linear SDE for which the usual Euler scheme diverges. 
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1 Introduction 

Let (f2, J 7 , P, Ft) be a complete probability space with a filtration and let a Wiener process 
{W t )t>o defined on this space. Consider the following stochastic differential equation, 

x t = x a + / a(x s )ds+ / b(x s )dW s , (1) 
Jo Jo 

where a, b : R+ xl->I are measurable functions and xq such that is J^-measurable 
and square integrable. Suppose that this problem has a unique strong solution x t . 

Assumption A Assume that there exist f(x,y), g(x,y) : M x 1 -> R such that 
f(x,x) = a(x), g(x,x) = b(x). Furthermore, assume that for any R > and any 
\xi\, \x 2 \, \y\\, 1 2/2 1 < R we have 

|/(xi,yi) - f(x 2 ,y 2 )\ + \g(xi,yi) -g(x2,V2)\ < C R (\xi - x 2 \ + I j/i - 2/2 1), 

for some Cr depending on R and on a, &, /, g. 

Let = to < t\ < ... < t n = T and set A = — . For any t £ [tk,tk+i] consider the 
following stochastic differential equation, 

Vt=yt k + f(Vs,yt k )ds+ g(y s ,y tk )dW s , t £ [t k ,t k+ i] (2) 
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where we assume that SDE (2), in each step, has a unique strong solution which we 
denote by y t and yo = xq. Note, that the coefficients in (2) are random. Let us call (2) 
a "semi-discrete" numerical scheme because, in general, we discretize only a part of our 
original SDE. In practice wc will choose /, g such that this numerical scheme will have 
a known explicit solution, for example, the "semi-discrete" SDE (2) may be linear. Wc 
can write the above numerical scheme more compactly, 

Vt=x Q + / f(y s ,y§)ds+ / g(y s ,y s )dW s , (3) 
Jo Jo 

where 

s = t k , s e [t k ,t k+ i]. 

The solution y t of (3) depends on A and we should use a notation like this, y^, but we 
don't do here for simplicity. Note that we use the second variable in /, g to denote the 
discretized part of the original SDE. 

The SDE (3), i.e. the "semi-discrete" SDE, is not an implicit numerical scheme, 
because in order to solve for y t we have to solve a stochastic differential equation and not 
an algebraic equation. In our setting, we can not reproduce the known implicit numerical 
schemes but we can reproduce for example the Euler scheme. So the usual Euler scheme 
belongs to our setting choosing f(x,y) = a(y) and g(x,y) = b(y). Another interesting 
way to choose /, g is f(x,y) = — \b (y)b(y) + a(y) + |6 (x)b(x) and g(x,y) — b(x) (see 
PP). This comes from the fact that the following SDE, 

x t = x + ( l-b (x s )b(x s )ds + J b{x s )dW Sl 
Jo * Jo 

has a known explicit solution (see [TJ, p. 117). We can use also other, more sophisticated, 
SDEs with known explicit solutions as these described in [7J. Thus, our method here is 
more general than Q] because we can arrive to other SDEs, like linear SDES, choosing 
suitable /, g. The main advantage of our method is that we produce always explicit 
numerical schemes in contrast to other interesting but implicit methods (see for example 
0, @). 

Another suitable choice, that we are going to use in our example in this paper, is 
f(x,y) = ^^-x and g(x,y) = ^jf^-x. Then the resulting "semi-discrete" SDE will be a 
linear stochastic differential equation with known explicit solution. 

Let us point out that our main result and setting remains true for the multidimen- 
sional case but in order to present our method as simply as possible we will remain in 
the one dimensional case. 

In the second section we will state and prove a convergence result and next we will 
give an application. For a general study of the numerical analysis of stochastic differential 
equations one can see [7j. 

2 Convergence of the semi discrete numerical scheme 

In this section we shall prove that our numerical scheme converges to the true solution. 
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Theorem 1 Assume that Assumption A holds and that (2) has a unique strong solution 
in each interval. Suppose that 

E( sup \x t \ p ) < A, E( sup \y t \ p ) < A 

0<t<T 0<t<T 

for some p > 2, A > independent of A. Then 

lim E sup \y(t) - x(t)\ 2 = 0. 

0<t<T 

Proof. Set p R = inf{t € [0,T] : \x t \ > R} and t r = inf{i £ [0,T] : \y t \ > R}. Let 
6 R = mm{T R ,p R }. 

At first, let us estimate the following probability 



\t r < T) = E 



L {r R <T}- 



RP 



< 



A 

~RP' 



Therefore we can prove that P(r r < T or p R < T) < Using Young inequality 
we obtain, for any 6 > 0, 

2 \ 2P +1 SA (p - 2)2A 



E ( sup \y t -x t \ < E sup \y t A9 R - x tA e R \ 

\0<t<T J \0<t<T 



pS?- 2 RP 



/(•,0. 



We shall estimate now the term \xt/\e R —yt/\9 R \ 2 as follows, using Assumption A for 



x tA9 R ~ VtA9 R I — 

tA9 R ft/\0 R 

(f(x s ,x s )- f(y s ,y$))ds+ / (g(x s ,x s ) - g(y s ,y s ))dW s 



t/\9 R ptA9 R 

2 1 \f(x s ,x s ) - f(y s ,y s )\ 2 ds + 2\ (g(x s , x s ) - g(y s ,y s ))dW s \ 

lo Jo 

r tA9 R r tA9 
|2 , |_ „, 2 



<C R (\x s -y s \ z + \x s -y s \ z )ds + 2\ (g(x s , x s ) - g(y s , y s ))dW s \ A 

Jo Jo 

Note that C R will be different from line to line. We can write now, 

Sup \x tA0R ~ VtA9 R \ 2 < Cr / (\x rA R - y T A9 R \ 2 + \XrA9 R - y^Jz\ 2 ) dr + 
0<t<s Jo 

ptA9 R 

2 sup | / (g{x s ,x s ) - g(y s ,y s ))dW s \ 2 . 

0<t<s Jo 

Taking expectations on both sides, using Doob's martingale inequality for the second 
term at the right hand side and Assumption A for g(-, ■) we arrive at, 



0<t<s 



E( sup \x tA 9 R - 2/tAe R | 2 ) < CrE / (\x rA e R - VrA9 R \ 2 + \x rA e R ~ V^0~^\ 2 ) dr ^ 



■f 

Jo 



Cr I E Q Sup \xiA0 R -yiAe R \ dr + C R j E|y rA9ji - J/^eji^' 



We shall estimate the term K\yt^g R — y^~ R \ 2 ■ We begin with, 

ptA6 R rtA0 R 

\ytA6 — VffjTr I 2 = I f(.Vs,ys)ds+ I g{y S7 y s )dW s \ 2 < 

JtA6 R Jtf\6 R 
tA6 R ftA6 R 

2 i ol / „. \JTJ/ |2 



2(/ f(y s ,y s )dsy + 2\ / g{y s ,y s )dW s \ 

JtA8 R JtA9 R 

Taking expectations, using Ito's isometry and the fact that \f(y s ,ys)\, \g(ys,ys)\ < Cr 
we have that, 

E\y tAS - yiKt) \ 2 <C R A. 
Collecting all together and using the Gronwall inequality we arrive at, 

2P +1 5A (p - 2)2A 



sup \y t -xtY) < C R A 
o<t<T ) P P 5~RP 

Now, given e > choose S > such that 2P+ p 5 A < s/3 and then choose i? such that 
(p-2)2A < £ ^3^ pinaiiy^ choose A small enough to get the desired result. 

□ 

3 Example 

Below we will propose an explicit numerical scheme for a super linear SDE as an appli- 
cation of Theorem 1. Consider the following SDE, 

xt = x - / x^ds + / bx s dW s , 
Jo Jo 

with xo £ IR+ and 6 g R. This SDE has a unique strong solution and we know that the 
usual Euler scheme diverges (see [3]). 

We propose the following "semi-discrete" numerical scheme which is explicit, 

t pt 

2-\ 



y t = x + / y s (-y s )ds+ / by s dW s . 
Jo Jo 

This semi discrete scheme is a linear SDE with unique strong solution, 

y t = xoe-^l+ b 4)ds+bW^ 

We need the moment bounds which we will prove in the next lemma. 
Lemma 1 Suppose that Xq > and Xq G K. Then there exists some A > such that 

E( sup \y t \ p ) < A, E( sup \x t \ p ) < A, E( sup -— ) < A, 

0<t<T 0<t<T 0<t<T \Xt\ 

for any p > 2. 
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Proof. Set r = xaxa{p,T} where p = mi{t e [0, T] : \x t \ > R} and t = mt{t 6 [0,T] : 
\y t \ > R} for some R > 0. 

Using Ito's formula on \ytr\ r \ p we obtain, 

\Vtr,r? = |X | P + f {p\Vsn-yi) + ^ \y s \ P )h(0,r)}( S ) dS + / bp|y»| P l[{(0.r)}(«)^ < 

|x | p + X) j \y s \n mr)} (s)ds + J o bp\y s \n mr)} (s)dW s 

Taking expectations we obtain, 

b 2 n(n-l) f l 
E|y t A,| p < W + ^ E|y sAr |^ S . 

Therefore, using Gronwall inequality, we arrive at, 

ny t Ar\ p < C(p)\x \', (4) 

with C(jp) independent of R > 0. But E|y tAr |P = E(|?/ tAr | I 'I {T .> t} ) + RPP(r < t). That 
means that P(t Ar < t) = P(r < t) —> as R —> oo so t Ar —> t in probability and noting 
that r increases as R increases we have that t A r — > t almost surely too, as R — > oo. 
Going back to (4) and using Fatou's lemma we obtain, 

< c( P )|x r, 

for any p > 2 and this is crucial to ensure that E J Q bp\y s \ p dW s = in the next step. 
Using Ito's formula again on |yt| p , taking supremum and then expectations and finally 
Doob's martingale inequality we arrive at, 

E( sup \y t \ p ) < A, 

0<t<T 

for some A > 0. To prove now that 

E( sup \x t \ p ) < A 

0<t<T 

we use Theorem 2.4.1 of [9] and obtain first the bound, 

E\x t \ p <A, 

for any p > 2 and then using Ito's formula on \x t \ p taking supremum and then expecta- 
tions and finally Doob's martingale inequality we obtain the desired bound. 
Finally, we shall prove that 

E( sup \) < A. 

0<t<T x t 

Set I = infjt > : — > ml. Using Ito's formula on (—^—) 2 we have, 
E( — ) 2 =E( — ) 2 +E / (2 + 3(— ) 2 )ds. 

XtAl %0 Jo x sf\l 
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Using Gronwall's inequality we can prove that E(^rL-) 2 < A with A independent of to. 
As before we deduce that E(^-) 2 < A and then that 

E( sup \) < A. 

0<t<T X t 

□ 

Now we can apply Theorem 1 to prove that our explicit numerical scheme converges 
in the mean square sense. Let us note that from the above moment bounds for the true 
solution we deduce that xt > a.s. Indeed, we have proved that 

E(^j-) = E(^-I {J>t} ) + m 2 P{l <t)<A. 

X tM X tAl 

Therefore, letting to — > oo we obtain that P(l < t) — > 0, noting that 

oo 

P(x t < 0) = P( Pi {x t < —}) = lim P({x t < —}) < lim P(l < t) = 0. 

' ' in n— »oo to m— >oc 

m— 1 

Thus, we need our numerical scheme to be positivity preserving. As we can see easily 
our scheme has this advantage. For this kind of problems there is also the tamed Euler 
scheme as the authors in [2] proposes. This kind of method behaves very well and in 
comparison with our "semi-discrete" method seems that the tamed Euler method is less 
expensive because we have in each step to compute an exponential. However, the tamed 
Euler scheme does not seem to preserve positivity, at least, for any A > 0. 

We shall give a simulation to show that the two numerical schemes are close. So, 
let the following SDE, 



x t = 1 + f ~x 3 s ds + I x s dW s , t G [0, 1] 
Jo Jo 



We apply the tamed Euler scheme and the semi discrete scheme to this problem for 
A = 1CP 4 and we plot the difference between these methods, i.e. if z t = yj amed — y^ eml 
we plot z t on [0, 1]. 



4 Conclusion 

In this paper we propose a new numerical method for solving stochastic differential 
equations. We apply our method to a super-linear SDE and compare with the tamed 
Euler method. We sec that our method preserves positivity of the true solution. Our 
method seems to behave very well when applied on super-linear problems. In order to 
manage other interesting problems a suitable transformation and then the application 
of this "semi-discrete" method maybe the answer. Let us point out that our main result 
remains true also for a multidimensional case. Our goal in the future is to apply our 
method to other stochastic differential equations arising in financial mathematics and 
to give a more detailed study of this method studying the rate of convergence, the case 
where the choice /, g are not locally Lipschitz etc. 



G 



Figure 1: Difference between the semi discrete scheme and tamed Euler scheme for 
x = l, A = 1CT 4 , 6 = 1, T = l. 
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Figure 2: Tamed Euler method does not preserve positivity, Xo = 1, A = 10 3 , b = 20, 
T = 1. 
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